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Experiments at the Large Hadron Collider (LHC) have measured multiplicity distributions in p+p and p+Pb 
collisions at a new domain of collision energy. Based on considering an energy-dependent broadening of the 
nucleon’s density distribution, charged hadron multiplicities are studied with the phenomenological saturation 
model and the evolution equation dependent saturation model. By assuming the saturation scale has a small 
dependence on the 3-dimensional root mean square (rms) radius at different energy, the theoretical results are in 
good agreement with the experimental data from CMS and ALICE collaboration. The predictive results in p+p 


collisions at \/s = 14 TeV of the LHC are also given. 
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I. INTRODUCTION 


Charged hadron multiplicities in relativistic heavy ion col- 
lisions are of great significance as its variables are very sen- 
sitive to the underlying mechanism involved in the nuclear 
collisions [1].These can provide a unique opportunity to test 
the predictions of quantum chromo-dynamics (QCD) and un- 
derstand the partonic structure of the colliding nuclei. Among 
the variables, the pseudo-rapidity distribution and negative bi- 
nominal distribution (NBD) are two important and interesting 
variable quantities to be studied on both experimental and the- 
oretical sides in recent years [2-7]. 

The pseudo-rapidity distribution of charged hadron is 
one quantity always used to investigate the particle produc- 
tion mechanism in high energy hadron-hadron and nucleus- 
nucleus collisions. Since pseudo-rapidity density is sensitive 
to the initial conditions of system and the hadronic final state 
interactions, the study of the charged hadron densities at mid- 
rapidity can provide the relevant information on the interplay 
between hard parton-parton scattering process and soft pro- 
cess. The NBD is another quantity which plays major role in 
describing multiplicity distribution of produced charged par- 
ticles. The NBD can describe well the multiplicity distribu- 
tion in almost all inelastic high energy collision processes ex- 
cept for the data at the highest available collider energies. In 
this paper, the pseudo-rapidity distribution and the NBD in 
p+p and p+Pb collisions will be studied in the framework of 
the Color Glass Condensate. 

At very high energies or small Bjorken-x, QCD predicts 
that high density gluon in a hadron wave function form a new 
state, the Color Glass Condensate (CGC). In this regime, the 
gluon density increases inside the hadron wave functions and 
becomes very large in comparison to all other parton species 
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(the valence quarks), and the sea quarks are suppressed by the 
coupling a, since they can be produced from the gluons by 
the splitting g > qq. 

Now, CGC has become an effective theory in describing 
the gluon content of a high energy hadron or nucleus in the 
saturation regime and many saturation models have been es- 
tablished [8-12]. These saturation models can be divided into 
two main kinds: the phenomenological saturation model and 
the evolution equation dependent saturation model. The phe- 
nomenological saturation models, such as the Golec-Biernat 
and Wusthoff (GBW) model [8], the Kharzeev, Levin, and 
Nardi (KLN) model [9], and the Kovchegov, Lu, and Reza- 
eian (KLR) model [10], are very simple and their dipole- 
proton scattering amplitude can be given by an analytic ex- 
pression. The evolution equation dependent saturation mod- 
els, such as the impact parameter dependent saturation model 
(IP-Sat) [11] and the running coupling Balitsky-Kovchegov 
(rcBK) model [12], are based on the DGLAP or the rcBK 
evolution equation and they are very useful at small Bjorken- 
x or high collision energy. Having studied charged particle 
multiplicity distribution with the GBW model [7], we will in- 
vestigate charged hadron pseudo-rapidity distribution and the 
NBD with the phenomenological KLN model and the evolu- 
tion equation dependent rcBK model, respectively. 


In order to give an accurate theoretical analysis, the nu- 
cleon’s density distribution in position space, which controls 
the unintegrated gluon distribution (UGD) through the satu- 
ration scale Q., must be considered. It should be pointed 
out that, due to gluon saturation, the width of the gluon dis- 
tribution inside a nucleon should grow with collision energy 
s/s [13]. This will lead to a broadening of the nucleon’s den- 
sity distribution in position space as ,/s increases. 

Thus, we consider an energy-dependent broadening of the 
Gaussian nucleon thickness function, and assume the satura- 
tion scale has a small dependence on the 3-dimensional root 
mean square (rms) radius at different collision energies [14]. 
The outline of this paper is as follows. The theoretical method 


is given in Sec. II and the results and discussion are given in 
Sec. I. 
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Il. METHOD 


For hadron-hadron collisions, the multiplicity per unit ra- 
pidity and per unit transverse area can be given by [15] 
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where N. = 3, £12 = (pz/\/8)exp(+y), b is the impact fac- 
tor and r is the transverse position of the m The running 


coupling constant &s(k?) = min GESTIN p0 5 


with 6o = 11 — Bnp = = 9 and A = Agcp = 0.2 GeV. For the 
unintegrated gluon distribution ¢, we will use the form given 
by the KLN model [9] and the rcBK model [12]. 

In the KLN model, ¢ is taken to be 


(zx, k?, b) = 
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where Cp = (N? —1)/(2N-) and « is a normalization fac- 
tor. The unintegrated gluon distribution depends on trans- 
verse position through the saturation scale 
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where Q3 = 2GeV?, A = 0.288 [8], and Tp,o is taken as 
1fm~?. For the nucleon thickness function, the Gaussian 
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and the kernel for the running term using Balitsky’s prescrip- 
tion reads 
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For the initial conditions, the GBW ansatz is used for the 

dipole scattering amplitude [8] 
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form is used 


e7% /(2B) 
T (bB = —~——_. 4 
pb) = 5 4 
where the proton width parameter B can be computed 
from [13] 


Tin( Vs) 2 
B = f 5 
(V8) = e Ee’ 6) 
and cin(y/5) is the inelastic scattering cross section. In this 
paper, we assume the gluon saturation scale has a small de- 
pendence on the 3-dimensional rms radius of the proton [14] 


2 
Q? (V8) = Q? o(/a0)( ys, (6) 


where ô = 0.8 [14] and the 3-dimensional rms radius ryms = 
\/(r2) = V3B. In Table 1 we collect a few representative 
values. 


TABLE 1. The 3-dimensional rms radius for various collision 
energies. The values for oj, at LHC energies were reported in 
Ref. [6, 16, 17] 

v/s (TeV) 0.9 2.36 5.02 7 14 
Cin (mb) 52 60 67 70.45 76.3 
VB (fm) 0.603 0.648 0.685 0.702 0.730 
Trms (fm) 1.044 1.192 1.186 1.216 1.264 


In the rcBK model, the unintegrated gluon distribution can 
be obtained from the dipole scattering amplitude via a Fourier 
transform 


2 
$(a, k) = / Se N (0,2) = I T Jo(rk)N (e,r), 
(7) 


where Jo is the spherical bessel function of the first kind. The 
dipole scattering amplitude in the rcBK evolution reads [12] 


N(ri, Y)N(ra, Y)I, (8) 


where y = 1 and Q2, is the initial saturation scale squared. 

Now let us take into account the negative binomial distri- 
bution. Negative binomial distribution is a general property 
of multi-particle production process regardless of type of col- 
liding particles, and the negative binomial probability distri- 
bution for obtaining n charged particles in the final state is 
given as follows: 


T(k+n) ny” nTn 
P= QED a 
(") = Fagrint 1) Me k Gi 
where the mean multiplicity ñ can be calculated in the frame- 
work of CGC by integrating y in Eq. (1). The quantity k, 
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(a) s'"=0.9 TeV 
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(b) s'"=2.36 TeV 


Fig. 1. Pseudo-rapidity distribution of charged hadrons in p+p collisions at v/s = 0.9 TeV (a), 2.36 TeV (b), 7 TeV (c), and 14 TeV (d). The 
solid and dashed curves are the result of the KLN model with and without considering the rms radius dependent saturation scale, respectively. 
The dotted curves are the result of the reBK model. The data is from CMS [2, 3]. 


which is the fluctuation parameter, can be estimated as a func- 
tion of the saturation scale [6] 


N?2-1 
k= WS —Q5 Us Vs)oin(Vs), (12) 


where «’ is a normalization factor. Here, in order to get an an- 
alytic expression between Qs,p and y, Qs,p is reconsidered by 
substituting x12 = (Qs/,/s)e*" into Eq. (3) as in Ref. [8]. 


II. RESULTS AND DISCUSSION 


To evaluate the pseudo-rapidity distribution, Eq. (1) should 
be rewritten using the transformation 


1. 1/cosh?n + m2/p2, + sinhn 


y(n) = 5m , (13) 
\/ cosh?°n + m3 /p2, — sinhn 


and the Jacobian can be correspondingly written as 


Ios Oy _ coshn (14) 
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where mo is the rest mass of particle, which corresponds to 
the order of the scale Agcp. With a y? analysis of the ex- 
perimental data [18], the factor « in Eq. (2) is equal to 0.51 
and 0.89 for the theory with and without considering the rms 
radius dependent saturation scale, respectively. In Fig. 1, the 
theoretical results for pseudo-rapidity distribution of charged 
hadrons in p+p collisions at \/s = 0.9 TeV (a), 2.36 TeV (b), 
7 TeV (c), and 14 TeV (d) are shown. The solid and dashed 
curves are the result of the KLN model with and without con- 
sidering the rms radius dependence of the saturation scale, 
respectively. The dotted curves are the result of the rcBK 
model. The experimental data comes from CMS [2, 3]. It is 
shown that the theoretical results of the KLN model consid- 
ering the rms radius dependence are in good agreement with 
the experimental data. For the effective value of Bjorken-x 
in the rcBK model is 1071? < x < 0.01, it is not valid at 
small \/s or large pseudo rapidity. Thus, only the results at 
vs = 7TeV fit well to the data with the rcBK model. The 
results shown in Fig. 1(d) are the predictive results for forth- 
coming LHC experiment at ys = 14TeV, and it is shown 
that the predictive results of the reBK model are almost the 
same as those of the KLN model. 


Figure 2 shows the negative binominal distribution of a 
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Fig. 2. Negative binominal distribution in p+p collisions at v/s = 0.9 TeV (a), 2.36 TeV (b), 7 TeV (c), and 14 TeV (d). The symbols and 


lines are the same as that in Fig. 1. The data is from ALICE [4]. 


charged hadron at |7| < 0.5. The symbols and lines in Fig. 2 
are the same as that in Fig. 1. The experimental data comes 
from ALICE [4]. The agreement is seen to be quite well for 
the KLN model considering the rms radius dependent satu- 
ration scale at \/s = 0.9 TeV and 2.36 TeV and for the reBK 
model at \/s = 7 TeV. Here, it should be noted that the agree- 
ment is seen to be not quite well for the theoretical results of 
the KLN model at ys = 7TeV even after considering the 
rms radius dependent saturation scale. Thus, we will give a 
systematic analysis of all 7 regions by considering the impact 
factor dependence of the mean multiplicity and the quantity 
k in the near future. 

Recently, the experimental data of the charged hadron mul- 
tiplicity in p+Pb collisions at \/s = 5.02 TeV are given by the 
ALICE collaboration [5]. In order to get a further test of the 
theory, the pseudo-rapidity distribution of charged hadrons in 
p+Pb collisions are investigated. For the nuclear density dis- 
tribution of Pb, we use the Woods-Saxon distribution [19] 


= PO 
p= 1 + exp|(r — R)/aļ’ aa 


where po corresponds to the nucleon density in the center of 
the nucleus, R is the nuclear radius and a is the “skin depth”. 
The theoretical results are shown in Fig. 3. The solid curve 
is the result of the KLN model and the dashed curve is the 
result that we gave in Ref. [7]. The dotted and the dash-dotted 
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Fig. 3. Pseudo-rapidity distribution of charged hadrons in mini- 
mum bias p+Pb collisions at ys = 5.02 TeV. The curves are the 
result of the KLN model (solid curve), the result that we gave in 
Ref. [7] (dashed curve), the results given in Ref. [20] (dotted curve) 
and Ref. [21] (dash-dotted curve). The data comes from ALICE [5]. 


curves are the results of DPMJET [20] and HIJING with the 
gluon shadowing parameter s, = 0.28 [21], respectively. It 
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is shown that the theoretical results of the KLN model fit well 
to the new experimental data. 


IV. CONCLUSION 


In summary, the pseudo-rapidity distribution and the NBD 
of a charged hadron in p+p and p+Pb collisions are studied 
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